Magnetic Properties of Ab initio Model for Iron-Based Superconductors LaFeAsO 
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By using variational Monte Carlo method, we examine an effective low-energy model for LaFeAsO 
derived from an ah initio downfolding scheme. We show that quantum and many-body fluctuations 
near a quantum critical point largely reduce the antiferromagnetic (AF) ordered moment and the 
model not only quantitatively reproduces the small ordered moment in LaFeAsO, but also explains 
the diverse dependence on LaFePO, BaFe2As2 and FeTe. We also find that LaFeAsO is under large 
orbital fluctuations, sandwiched by the AF Mott insulator and weakly correlated metals. The orbital 
fluctuations and Dirac-cone dispersion hold keys for the diverse magnetic properties. 
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Discovery of iron-based superconductors has renewed 
interest on mechanisms of high- Tc superconductivity [l[ . 
Roles of electron correlations in the iron-based fami- 
lies and in the pairing mechanisms are under strong de- 
bates An aspect common in these families is the 
nearby antiferromagnetic (AF) phases similarly to the 
cuprate high- Tc superconductors. In the iron-based fam- 
ilies, however, the AF orders are found in metals and 
show a variety of ordered moment ranging from ^ 2 jUb 
with a bicollinear order for FeTe to 0.3-0.6 ub with an 
AF stripe (AFS) order for LaFeAsO [13, Eli. Since ah 
initio density-functional calculations usually underesti- 
mate AF moments, it is unconventional that the mea- 
sured moment for LaFeAsO is substantially smaller than 
the density-functional estimate of 1.2-2.6 ^b [IMl. It 
strongly suggests the necessity of considering quantum 
fluctuations beyond the mean-field level. 

A difficulty in the iron-based superconductors comes 
from an entangled band structure composed of five Fe-3d 
bands near the Fermi level [l^ ■ Recently proposed three- 
stage scheme, consisting of the global band structure by 
the conventional density functional calculations, ab ini- 
tio downfolding scheme to derive low-energy models, and 
solving the resultant realistic models, has opened a way 
of analyzing such a real complexity of materials 16|, |17 1 . 
This scheme has already been applied to the iron-based 
families and effective low-energy models have been de- 
rived [l^, [3]. The next step of solving the effective 
models by reliable low-energy solvers is so far mostly con- 
fined to that by the dynamical mean-field approximation 
(DMFA) SHU, where spatial correlation effects are 
hardly analyzed in the present multi-orbital systems. 

To understand the correlation effects and the uncon- 
ventional magnetism described above, interplays of or- 
bitals and spins have to be elucidated by considering spa- 
tial fluctuations beyond DMFA. For this purpose, many- 
variable variational Monte Carlo (VMC) method [20|, 
combined with quantum-number projection (2l| . offers 
a suitable and accurate low-energy solver in clarifying 
dynamically and spatially fluctuating phenomena [20|. 



In this letter, we apply unrestricted Hartree-Fock 
(UHF) and VMC methods to solve the downfolded ef- 
fective low-energy model of LaFeAsO (3]. We clarify 
the microscopic origin of the small AF ordered moment. 
In addition, we find that LaFeAsO is located in a region 
of large orbital fluctuations, sandwiched by the AF Mott 
insulators and weakly correlated metals. 

Our low-energy model derived from the ab initio 
scheme [l3| is defined for ten-fold degenerate Fe-Sd or- 
bitals in a unit cell containing two Fe atoms in the form 
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where {^ur) ^ creation (annihilation) operator of 
an electron with smn a in the nth maximally localized 
Wannier orbitals [2^ centered on Fe atoms in the unitcell 
at R. tmRnR' contains single-particle levels and transfer 
integrals. f/mRjiR' and JmRnR' are screened Coulomb 
and exchange interactions, respectively. Offsite interac- 
tions were dropped since those are more than four times 
smaller than the onsite parameters. Details of the pa- 
rameter derivation are given in Ref. [l8l | . 

Here, we introduce a measure of the geometrical frus- 
tration defined as a ratio of the diagonal next-nearest 
transfer [t') to the nearest-neighbor transfer (i), i.e., t'/t. 
For yz/zx and — orbitals, this value is roughly 
1.0, while for xy and orbitals t' /t - 0.1 Thus, 
yz/zx and x'^ — y^ orbitals are categorized to strongly 
frustrated orbitals whereas xy and orbitals arc unfrus- 
trated. This distinction in the degree of frustration is a 
characteristic feature of this model and controls magnetic 
properties. Another point to be stressed is an apprecia- 
ble orbital dependence of the Coulomb interaction; the 
smallest intra-orbital Coulomb interaction is U^'^-y^ = 
2.20 eV and the largest one is Uxy = 3.31 eV. Such dif- 
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FIG. 1: (color online) Left panel: Ground-state UHF energy 
per site of AF and paramagnetic metal (PM) as functions of 
interaction ratio A. Right panel: Four AF patterns. Energy 
is measured from that of AFG. The ground state is a PM for 
< A < 0.5, AFS for 0.5 < A < 1.1, and AFG for A > 1.2. 
Around A ~ 1.1, AFH is the ground state. For A < 0.5, all 
the solutions converge to PM. 



ferences bring about different roles for different orbitals, 
namely an orbital differentiation by lifting the degener- 
acy, in determining magnetic and charge structures. 

In the multi-band model, effective one-body potential 
li is given by fiv = fJ-u + J2„' U^'rii,' + where 
/i^ is a one-body potential for orbital v and is its 
occupation. The Hartree term described as the second 
and third terms in is already included in the LDA 
calculations. Therefore, we need to eliminate this Hartree 
contribution to exclude the double counting. To this end, 
we correct ^i^, so as to satisfy /i^ = where 
is the one-body potential obtained from LDA. For multi- 
band systems this correction is necessary. On the other 
hand, we ignore the double counting from the exchange 
correlation energy by LDA, because it is small [24 1. 

In order to monitor the Coulomb interactions U as 
well as the exchange interactions J, we introduce the 
interaction ratio A to scale all the matrix elements of 
U and J uniformly. The original ab initio values (non- 
interacting case) correspond to A = 1 (A = 0). 

To gain a rough insight into the effective Hamiltonian 
in Eq. ([T|), we show results of UHF calculations in Fig. [T] 
Although the UHF overestimates stability of the ordered 
state, it is useful to understand the global structure of the 
model. Wc consider four different AF solutions as shown 
in the right panel of Fig. [1] AFS, G-type AF state (AFG), 
and bicoUinear AF state (AFB) as the candidates for 
the ground states of the iron pnictidcs. In addition, we 
examine the stability of the half-collinear AF state (AFH) 
known to be realized in a frustrated Hubbard model [2lj . 

In the left panel of Fig. [U we plot their energies as 
functions of A. Here, is fixed at the value at the 
original interaction parameter ( A = 1). With increasing 
A, paramagnetic metal undergoes an AF phase transition 
around A = 0.5 into the AFS state. The AFS phase 



continues to be stable through A = 1. For A > 1.2, 
AFG becomes the ground state. Sandwiched by AFS 
and AFG, AFH is stabilized around A ~ 1.1, implying a 
relevant frustration effect in this region. The AFG phase 
is always the ground state when we ignore the double 
counting correction of the Hartree terms. 

Orbital occupations rii, are monitored as functions of 
A in Fig. [2] (a). From this plot, we can understand why 
the AFS (AFG) phase is stable for A < 1.2 (A > 1.2): 
For A ^ 1, all the orbitals contribute to the magnetic 
moment; each orbital occupation is close to half filling 
{rii, ^1). In this region, the number of the frustrated 
orbitals {yz/ zx, — if') is larger than that of the un- 
frustrated orbitals [xy, z'^). Hence, the former frustrated 
orbitals dominate the magnetic structure and AFS be- 
comes stable. In contrast, for A > 1.2, nj;2_y2 increases 
rapidly away from half filling and loses magnetic activity. 
Such a disruption in the subtle balance of the frustrated 
and unfrustrated orbitals leads the latter orbitals, having 
the stronger Coulomb interaction {Uxy ~ 3.31 eV, 11^2 = 
3.27 eV) than the former {Uy;, = U^a: = 2.77 eV), to dom- 
inate the magnetic structure; AFG is realized for large A. 
The result observed here is interpreted as electron differ- 
entiation in orbitals; the appreciable orbital dependence 
of the frustration and interaction parameters generates 
emergent distinction in roles of each orbital. 

Figure [2] (b) plots A dependence of the mag- 
netic ordered moment defined by m{q)'^ = 
' E^.AS^SJ)e"l^'''-''^^ for AT.-site (iV,/2 unit 



cell) system with the periodic boundary condition and q 
is set to Qpeak = (0,7r). In this definition, the saturated 
magnetic moment for the classical Neel state is given by 
m = 4 ^B- We also superpose density of states (DOS) 
at the Fermi level [p{E-p)] in the same plot. The AF 
phase transition occurs around A = 0.5 and then the 
metal-insulator transition occurs around A = 0.8. The 
AF metal thus exists for 0.5 < A < 0.8. Notice that 
the metal survives even for an unexpectedly large AF 
moment, m ^ 3.5 close to the saturated value. 

Although real LaFeAsO corresponds to A ~ 0.5 on 
the UHF level, it is illuminating to show UHF band 
structures at a large A = 0.7 near the metal-insulator 
transition to understand the origin of the robust metal 
(sec Fig. [2] (c)). Near the Fermi level, two Dirac nodes 
(marked as dots) appear along the M-F line 25, 2^. Un- 
til their pair annihilation, the Dirac nodes maintain the 
metallic band and the metal becomes stable even when 
a large AF ordered moment grows at a large A (~ 1). 

In Fig. [2] (d), we show partial DOS for A = 1.0 corre- 
sponding to the insulating region. Just below the Fermi 
level, the — y^ DOS is much smaller than those of the 
other orbitals, indicating a relatively large gap opening 
(^ 2 eV) for this orbital. This is consistent with the 
fact that, around A ~ 1, nx2_y2 is pinned to half filling 
while the other orbital occupations change continuously 
[sec Fig. [2] (a)]. It has similarity to the orbital selective 
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FIG. 2: (color online) UHF results for our ab initio low-energy 
model: (a) Orbital occupations as functions of A. Near A = 1, 
occupation of the x'^ — orbital is pinned at 1. (b) Mag- 
netic ordered moment m{qpca,k) and DOS at the Fermi level 
p{Ef)- (c) Band structures for A = 0.7 along the high- 
symmetry points, where F = (0,0), M = (ti",?!"), Y = (0, tt), 
and M = (— 7r,7r). Energy zero is the Fermi level. Dirac nodes 
(dots)appear along the M-F line, (d) Partial DOS at A = 1.0. 



Mott transitions [27[ in the sense that the charge gap 
depends largely on the orbitals, although here, all the 
orbitals have nonzero charge gap. 

To examine how quantum fluctuations beyond the 
mean-field approximation affect the magnetic properties, 
we performed VMC calculations. Our variational wave 
function |2(| is defined as 



(2) 



where Vg is the Gutz wilier factor. The spin (momen- 
tum) quantum-number projection C^^'^ (C^^^) restores 
the SU(2) spin-rotational (translational) symmetry and 
generates a state with the correct total spin 5=0 (to- 
tal momentum K=0). The one-body part Impair) is the 
generalized pairing wave function defined as li/ipair) ~ 
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FIG. 3: (color online) Benchmark of VMC on two-orbital 
model for A's =4x2. (a) A dependence of total energy 
obtained by UHF and VMC, measured from the exact result 
-Eed- (b) A dependence of magnetic ordered moment. 



VMC. In LaFeAsO, because doubly degenerate yzj zx or- 
bitals have the largest weight on the Fermi surface, we 
consider a two-orbital model for these orbitals with the 
same interactions as the original five-orbital models. This 
two-orbital model is strongly correlated and highly frus- 
trated; i.e., t' jt - 1.0 and IJ jt - 10. 

We compare in Fig. [3] (a) UHF, VMC, and exact- 
diagonalization (ED) energies of an iVs=4x 2 system. The 
ground-state energy of VMC becomes much lower than 
the UHF energy and very close to the ED one. Fig- 
ure[3](b) illustrates that quantum and many-body fluctu- 
ations drastically decrease the ordered moment m(qpcak) 
and VMC reproduces satisfactorily the ED result from 
the weak to strong coupling regimes. Clearly, our VMC 
works well even for serious regions of frustrations and 
electron correlations in the multi-orbital system. 

We give in Fig. |4] the VMC results for the low-energy 
effective Hamiltonian for LaFeAsO. The VMC ground 
state is found to be AFS for A < 1. System size de- 
pendence of m((3'poak) for various choices of A is shown 
in Fig. U (a). We extrapolated the data by l/v'TVs, in 
which finite size effects arc assumed to arise from linear 
spin-wave type excitations in the AF ordered states. The 
extrapolation gives to > only for Ac > 0.75, while, for 
A < 0.75, the ground state is paramagnetic. 

Here, there exist two origins of the overestimate of the 
interaction parameters derived in Ref. [T^]. First, Ref.[l8l| 
ignored the screening effects from La-/ electrons. This 
screening was shown to lead to the reduction of the aver- 
age of orbital diagonal Coulomb interactions from C/avp to 



|0) with Uj being the variational 0.9J7ave The second is the interlayer screening effect, 



parameters. In this study, we allow /.y to have 2x2 sub- 
lattice structure or equivalently we have 2 x 2 x 5^ x A'g 
parameters. All the variational parameters are simulta- 
neously optimized by using the stochastic reconfiguration 
method [20|, |28[. The variational function |?/;) in Eq. ([2]) 



can describe paramagnetic metals and the AFC and AFS 
phases, as well as superconducting phases. 

Before applying the VMC method to our model, we 
show a benchmark demonstrating high accuracies of our 



which has to be appended when we employ purely two- 
dimensional models as in the present study. The screen- 
ing by electrons on neighboring metallic layers reduces 
C/avo to 0.85C/avc for LaFeAsO [2§|. On closer inspec- 
tion, [/avc decreases to ^ 0.75J7avo in total. Therefore, to 
compare with the experimental results, we should simu- 
late correlation effects by A = (C/avc/C^avc)A instead of A, 
where C/ave is the corrected interaction J7avo ~ 0.75J7avc- 
The A dependence of to ((/peak) is shown in Fig. 0] (b). 



4 



(a) 
4 



^3 







(b) 



□ x= 


= 1.0 <> x= 


=0.76 




-0.9 


=0.75 


-Q X= 


=0.8 A X- 


=0.7 


X- 


=0.775 O X- 


=0.6 






♦ 








1 


1 1 





0.5 



^1 



BaFe2As2 



LaFeAsO 



LaFePO 



i 







0.2 0.4 

l/iNs) "2 



<^-^ 



0.8 




FIG. 4: (color online) Magnetic ordered moment 7Ti(gpoak) 
calculated by VMC for the model of LaFeAsO. (a) Size de- 
pendence of m((jrpcak) for several A. (b) A dependence of 
'Ti(gpcak) in the thermodynamic limit (open circles). Exper- 
imentally observed materials dependence at corresponding A 
is also shown by crosses. Quantum critical point of the AF 
transition appears at slightly below A = 0.75(A ~ 1). 



It indicates that LaFeAsO is close to the AF quantum 
critical point consistently with experiments [30|. Al- 
though details of the model parameters are not fully im- 
plemented, the ab initio parameters [31 tell us that other 
compounds are also simulated roughly by the correspond- 
ing scaled ratio of A specified from the ratio of Uavc/t 
with the averaged transfer t in Ref.[l^. Estimated or- 
dered moments for LaFePO, BaFs2As2 and FeTe are also 
all consistent with the experiments 13, lH [30l - l33l | . The 
relevance of quantum criticality is a recent subject of de- 
bates (33 - [36| . The present work supports its relevance. 

Finally, we describe why quantum fluctuations largely 
reduce the ordered moment. Around A = 1, each orbital 
occupation is at incommensurate filling, slightly away 
from half filling. This is analogous to the doped Mott 
insulators when we consider each orbital occupation, al- 
though the total filling is always commensurate. As a re- 
sult of this incommensurate filling, quantum fluctuations 
in orbitals substantially destroy the AF ordered moment. 

In summary, we have studied the magnetic properties 
of the ab initio low-energy model for LaFeAsO. State-of- 
the-art VMC calculations show that the realistic param- 
eter stabilizes a stripe-type AF metal close to a quantum 
critical point in agreement with the small ordered mo- 
ment observed in the experiment. Furthermore, a sharp 
A dependence of the ordered moment quantitatively ex- 
plains the ordered moment of LaFePO, BaFe2As2 and 
FeTc as well on unified grounds. Large fluctuations and 
differentiations in orbitals coexisting with a robust metal- 
lic dispersion arising from the Dirac cones also hold the 
key to understanding the low-energy physics. 

The authors gratefully thank Daisukc Tahara for the 
use of his VMC code and useful comments. 
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